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ABSTRACT 

In this work we test the assertion that the scatter in the mass of black holes which 
drive quasars should be luminosity dependent with less scatter in more luminous 
objects. To this end we measure the width of the Mgii A2799 line in quasar spectra 
from the SDSS, 2QZ and 2SLAQ surveys and, by invoking an unnormalised virial mass 
estimator, relate the scatter in line width to the scatter in mass in the underlying black 
hole population. We find conclusive evidence for a trend such that there is less scatter 
in line width, and hence black hole mass, in more luminous objects. 

However, the most luminous objects in our sample show such a low degree of 
scatter in line width that, when combined with measures for the intrinsic scatter in the 
radius-luminosity relation for the broad-line region in active galaxies, an inconsistency 
arises in the virial technique for estimating black hole masses. This analysis implies 
that, at least for the most luminous quasars, either there is little-to-no intrinsic scatter 
in the radius-luminosity relation or the Mg ii broad emission line region is not totally 
dominated by virial velocities. 

Finally we exploit the measured scatter in line widths to constrain models for 
the velocity field of the broad-line region. We show that the lack of scatter in broad 
line-widths for luminous quasars is inconsistent with a pure planar/disk-like geometry 
for the broad-line region. In the case of a broad-line region with purely polar flows 
the opening angle to luminous quasars must be less than ~ 55°. We then explore the 
effects of adding a random or spherically symmetric component to the velocities of 
gas clouds in the broad-line region. Assuming an opening angle to quasars of 45° a 
planar field can be made consistent with our results if '^ 40 — 50 % of the velocities 
are randomly distributed. 

Key words: galaxies: evolution - quasars: general - quasars: emission lines - cos- 
mology: observations 



1 INTRODUCTION 
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It is now apparent that the majority of massive galaxies 
harbour a super-massive black hole (SMBH) at their cen- 
tre. Dynamical studies of the sphere of influence of these 
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SMBHs have been successful in determining the mass for 
some tens of systems. In these cases the mass of the central 
SMBH has been observed to correlate strongly with proper- 
ties of their host spheroid such as luminosity (Magorrian et 
al. 1998), velocity dispersion {o\ Ferrarese & Merritt 2000, 
Gebhardt et al. 2000) and concentration/Sersic index (Gra- 
ham et al. 2001). More recently, Hopkins et al. (2007a) have 
demonstrated that these relationships can all be regarded 
as various projections of a "black hole fundamental plane" 
(BHFP), relating black hole mass to the potential well of 
the galaxy. 

The existence of these correlations and, in particular, 
the similarity of the BHFP to the fundamental plane of ellip- 
tical galaxies points to an intimate link between the growth 
of SMBHs and galaxy evolution. Moreover, the BHFP is con- 
sistent with the notion that SMBHs grew in a self-regulated 
manner, from gas gravitationally confined in the galaxy cen- 
tre, which was eventually expelled by feedback processes 
(Hopkins et al. 2007b). Recent hydrodynamical simulations 
that incorporate radiative cooling, star formation, black hole 
growth and feedback from both supernovae and nuclear ac- 
tivity have shown that mergers between gas-rich galaxies 
of comparable mass provide a mechanism for concentrating 
gas in galaxy centres through tidal effects (e.g. Barnes & 
Hernquist 1991, 1996), fueling the growth of SMBHs (e.g. 
Di Matteo et al. 2005), reproducing the correlations (e.g. 
Robertson et al. 2006), and explaining observed properties 
of quasars as a phase of evolution prior to the termination 
of black hole accretion (e.g. Hopkins et al. 2005a,b,c). 

Further evidence for a merger-driven origin of quasars 
is provided by comparing the evolution of the abundance 
of luminous quasars with that of the cosmic star forma- 
tion rate (SFR). From an empirical determination of the 
bolometric quasar luminosity function, Hopkins, Richards 
& Hernquist (2007) infer that the quasar luminosity den- 
sity peaks at z = 2.15 ± 0.05 and rapidly declines towards 
higher redshifts. Some observational estimates suggest a sim- 
ilar behaviour for the evolution of the cosmic SFR (e.g. Fan 
et al. 2001; Hopkins & Beacom 2006). However, at high red- 
shifts, incompleteness, cosmic variance of the surveys, and 
uncertain corrections owing to dust extinction complicate 
this analysis. From a measurement of the opacity of the 
Lyman- alpha forest (Faucher-Giguere et al. 2007a), Faucher- 
Giguere et al. (2007b) have shown that the optical depth at 
« = 4 is incompatible with a steep decline in the cosmic 
SFR at z > 3, as suggested by e.g. the results of Hopkins 
& Beacom (2006), but in accord with theoretical model- 
ing (e.g. Springel & Hernquist 2003; Hernquist & Springel 
2003) which predicts that the SFR should peak at z > 4. 
The implied offset between the peak in the cosmic SFR and 
the quasar luminosity density indicates that the growth of 
SMBHs is not directly tied to star formation or gas density, 
but is related to a secondary process. By employing esti- 
mates of the rate of halo mergers, Hopkins et al. (2007c) 
have shown that the quasar luminosity function can be re- 
produced if SMBH growth occurs primarily in major merg- 
ers of gas- rich galaxies. 

Motivated by these various lines of argument, Hop- 
kins et al. (2006a, b, 2007a,b) have developed a model 
for galaxy evolution in which mergers, starbursts, quasars, 
SMBH growth, and the formation of ellipticals are connected 
through an evolutionary sequence. As part of this work, Hop- 



kins et al. (2005a) used simulations of merging galaxies to 
quantify the phases of evolution associated with quasar ac- 
tivity and showed, in particular, that quasar lifetimes de- 
pend not only on the instantaneous luminosity of a quasar, 
but also its 'peak' luminosity. Convolving these lifetimes 
with estimates of the merger rates of galaxies, Hopkins et al. 
(2005b, c; 2006a; 2007a) were able to reproduce the observed 
optical and X-ray luminosity functions of quasars. 

In this interpretation of the quasar luminosity function, 
the brighter (> L*) objects are all massive black holes, ac- 
creting near the Eddington rate towards the end of their 
growth phase. Less luminous quasars can be either low-mass 
systems accreting rapidly and undergoing significant growth, 
or larger black holes accreting at comparatively lower rates. 

This model implies that the range of accretion rates in 
quasars should be luminosity dependent. That is, the range 
accretion rates (and by extension SMBH masses) at a given 
luminosity should be larger for lower luminosity objects, and 
decrease as the luminosity increases. 

In this work we aim to derive the dispersion of the 
SMBH mass of quasars as a function of luminosity. We mea- 
sure the Mgll broad emission-line width in spectra from 
three large spectroscopic surveys of Type I active galactic 
nuclei (AGN). Assuming there exists a virial relation for 
calculating SMBH masses from the Mgll line width we re- 
late the measured distribution of emission line widths to the 
SMBH mass distribution. 



2 VIRIAL SMBH MASS ESTIMATORS 

By far the most accurate, robust and believable method for 
measuring the mass of SMBHs is to perform dynamical stud- 
ies of stars, gas or masers in the potential of the SMBH (e.g. 
Herrnstein et al. 1999; de Francesco et al. 2006; Gebhardt 
et al. 2003). Rarely, however, is this viable for type I AGN 
which for the most part are too distant to resolve the sphere 
of infiuence of the SMBH, and have nuclei so bright that it 
is difficult if not impossible to obtain detailed photometric 
or spectroscopic information of the environment at redshifts 
- 1. 

The most direct method for calculating SMBH masses 
of type I AGN is via reverberation mapping of the broad-line 
region (BLR) and the virial theorem (Peterson et al. 2004). 
In this case the size of the BLR (rBLa) is given by the time 
delay between continuum and emission line variations and 
the velocity dispersion of the BLR (Vblr) is measured as 
the width of the emission line in the variable spectrum. The 
mass of the SMBH is then estimated as 



Mb 



/ 



G 



(1) 



The factor / defines what is not known about the BLR; 
its geometry, velocity field and orientation. The value of / 
is order unity and authors take various approaches assign- 
ing it a value. One can find / = 3/4 applicable to random 
orbits (Peterson & Wandel 1999), / = 1/(4 sin^ 61) for a disk 
inclined at an angle 6 to the observer (McLure & Dunlop 
2001), or even / = 1 for simplicity (McLure & Jarvis 2002). 
Onken et al. (2004) measured / = 1.4 by comparing the 
M — a relation in a group of reverberation mapped AGN in 
local spheroids. 
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This last value is somewhat higher than expected com- 
pared to the simple theoretical values. However, the inter- 
pretation of this is unclear and the discrepancy may be due 
to selection bias (Lauer et al. 2007) and/or cosmological evo- 
lution in the M — a relation (Woo et al. 2006) both of which 
could bias the measured value of / high. 

Due primarily to time constraints, the number of rever- 
beration mapped systems is in the tens and the luminosity 
range which these measurements span is not huge (Kaspi 
2007). However, one key result which has come out of rever- 
beration mapping is the radius-luminosity relation (Kaspi et 
al. 2005). While there is significant scatter in this relation, 
including significant intrinsic scatter of up to 40 %, it does 
offer a simple single epoch method for estimating the radius 
of the BLR in AGN. 

If one takes the instantaneous FWHM of a broad- 
emission line as a measure of the virial velocity in the BLR, 
then the r — L relation provides a method for estimating the 
unsealed SMBH mass with single epoch observations. This 
technique has become known as the virial method. To date 
virial relations have been calibrated for the H/3 (Kaspi et 
al. 2000; Vestergaard & Peterson 2006), Mgll (McLure & 
Jarvis 2002; McLure & Dunlop 2004; KoUmeier et al. 2006; 
Salviander et al. 2007; McGill et al. 2008) and Civ lines 
(Vestergaard 2002; Vestergaard & Peterson 2006). In each 
case these calibrations result in expressions of the form 



Mbh = A{\LxfFWHM^ 



(2) 



Where FWHM is the full width at half maximum of the 
spectral line in question and XL\ is the monochromatic lu- 
minosity of the continuum near that line. A is a normal- 
isation constant and the exponent a gives the luminosity 
dependence of the r — L relation. 



2.1 Problems with virial mass estimators 

When using virial SMBH masses it is prudent to state some 
caveats. Firstly and most obviously the virial relations are 
statistical in nature. Hence while they may be accurate when 
averaged over a large number of systems individual measure- 
ments should be viewed with caution. 

Secondly a number of studies have investigated the dif- 
fering emission regions for high and low ionisation lines (e.g. 
Richards et al. 2002; Baskin & Laor 2004; Elvis 2004 which 
raises questions as to whether all broad lines can be used as 
virial mass indicators. This is, however, still a subject open 
to debate. Vestergaard & Peterson (2006) reanalysed Baskin 
& Laor's data and showed that the apparent discrepancies 
between H/3 and C iv virial mass estimates are nullified by 
applying more appropriate selection criteria. 

Thirdly it is important to note that the calibrations of 
the Mgll and C iv virial relations are not direct. These lines 
have not been sufficiently studied in reverberation mapping 
campaigns to calibrate any potential r — L relation for them. 
Instead the virial relations for these lines are normalised 
through comparisons with SMBH masses calculated from 
the H/3 line. 

Finally, there is simply a dearth of solid information 
constraining the velocity field and geometry of the BLR. Re- 
verberation mapping as a technique has the potential to de- 
scribe the BLR in detail given sufficient quality data (Welsh 



& Home 1991). However, in reality this sort of idealised pre- 
cision is unlikely and to date these types of results have not 
materialised. Theoretical models for BLRs exist (e.g. Em- 
mering, Blandford & Shlosman 1992; Konigl & Kartje 1994; 
Murray & Chiang 1996) but the lack of strong observational 
constraints makes a proper comparison between these dif- 
ficult. However, almost all of these models imply a strong 
virial component to the velocity field of the BLR. 

There is good evidence that the virial method, while 
imprecise, can on average give accurate black hole masses. 
Comparisons between reverberation masses and bulge veloc- 
ity dispersion show an M — cr relation analogous to that ob- 
served in nearby quiescent galaxies (Onken et al. 2004). And 
comparisons between virial and reverberation masses also 
show agreement (Vestergaard & Peterson 2006; McLure & 
Jarvis 2002). But to obtain strong results from virial masses 
large samples are required to beat down random errors in 
these estimates. 

In this work we assume a virial relation for the Mgll 
line as given by equation 2. We do not, however, assume 
a parametrisation for this relation since it is not required. 
We take a large number of QSO spectra from several spec- 
troscopic surveys and bin them by luminosity and redshift. 
We invoke the r — L relation and assume objects in a lu- 
minosity bin will have the same tblr, and the scatter in 
SMBH masses within that bin is simply given by the scatter 
in broad-line widths. Furthermore, looking at the scatter in 
log space for sufficiently small luminosity bins we can ignore 
the L term in equation 2, and the coefficient and 



Disp(log(MBH)) = 2 X Disp(log(FW^HMMgii)) 



(3) 



Where Disp() denotes the dispersion in the given variable. 
Note that we have ignored all extraneous sources of scatter 
in this equation, some of which will be significant to our 
calculations. We discuss these in more detail in sections 6.3 
and 6.4. 



3 DATA 

We take as our sample all of the quasar spectra from the 
Sloan Digital Sky Survey (SDSS; York et al. 2000) data re- 
lease five (DR5; Adelman-McCarthy et al. 2007) as compiled 
by Schneider et al. (2007) . To increase our luminosity range 
we also take all QSO spectra from the 2df QSO Redshift 
survey (2QZ; Croom et al. 2004) and 2dF SDSS LRG And 
QSO survey (2SLAQ; Richards et al. 2005; Croom et al. in 
prep.) as well. Table 1 shows a brief summary of the number 
of objects and magnitude limits in each sample. 

We thus have spectra from three different surveys taken 
with two different instruments. We give here a brief descrip- 
tion of the spectra. 



3.1 SDSS spectra 

Details of the Sloan telescope and spectrograph are given in 
Gunn et al. (2006) and Stoughton et al. (2002). The spectra 
have a logarithmic wavelength scale translating to a disper- 
sion of ~ 1 — 2Apix~^ and a resolution A/ A A ~ 1800 in 
the wavelength range 3800 — 9200 A. Objects are observed 
initially for 2700 sec. Then are reobserved in 900 sec blocks 
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Table 1. Summary of the surveys from which we obtained spectra. Successive surveys have fewer spectra but go deeper, increasing 
our luminosity range at a given redshift. Note that the magnitude limits quoted for the SDSS QSO survey are those for the primary 
QSO survey. The high redshift sample goes deeper, and included are sources observed under different selection criteria and also QSOs 
identified as part of other surveys. 



Survey 



No. of Objects 



Mag. Limits 



Resolution Dispersion 



SDSS {DR5) 77,429 19 > i > 15 ~165km/s ~ l.sApix" 

2QZ 23,338 20.85 > 6,j > 18.25 ~465km/s ~ 4.3Apix- 

2SLAQ 8,492 21.85 > g > 18.00 ~465km/s ~4.3Apix- 



S/N 



^ 13pix-i 
5.5 pix^-*^ 
5.5 pix^-*^ 



until the median S/N is greater than -^ 4pix^^resulting in 
a S/N distribution with a mean at ~ 13pix~^. 

The spectra are extracted and reduced with the SPEC- 
tro2d pipeline and automatically classified with SPEC- 
troId (Stoughton et al. 2002). However, in creating the 
Sloan QSO sample used in this paper, Schneider et al. (2007) 
visually inspect all of the candidate spectra to determine 
their classification. 



3.2 2dF spectra 

Both 2QZ and 2SLAQ spectra were taken with the 2 degree 
Field (2dF) instrument on the Anglo-Australian Telescope 
with the 300B grating (Lewis et al. 2002). Spectra have a dis- 
persion of 4.3 A pix^^ and a resolution of ~ 9 A in the wave- 
length range 3700 — 7900 A. 2QZ observations were between 
3300 and 3600 sec compared with 14400 sec for 2SLAQ. The 
increase in exposure time for the fainter 2SLAQ sample re- 
sults in S/N distributions that are almost indistinguishable. 
Both peak at ~ 5.5 pix^^ for positive QSO IDs. 2dF spectra 
are extracted and manually classified during the observing 
run with the 2dFDR pipeline (Bailey & Glazebrook 1999) 
and AUTOZ redshifting code (Groom et al. 2001). 

The main difference between reduced Sloan and 2dF 
spectra is the lack of flux calibration for 2dF sources. An 
average flux calibration for the 300B grating has been cal- 
culated by Lewis et al. (2002) as part of the 2dF Galaxy 
Redshift Survey and in our analysis we do apply this correc- 
tion. However, a quick examination of 2dF spectra shows the 
inadequacy of this median calibration to correct for the vari- 
ations in response between differing spectra as illustrated in 
Fig.l. 



4 ANALYSIS 

Our goal is to measure the dispersion in the widths of broad 
Mgll lines in a large sample of QSOs. For this we require 
an automated routine to measure the width of Mgll in a 
consistent way across our sample and a robust method for 
analysing the results. In this section we discuss the line fit- 
ting code used to measure the line widths, and our analysis 
of the output, but first a quick note on absolute magnitudes. 



4.1 Magnitude calculations 

Magnitudes are given in the 6j-band in the 2QZ catalogue, 
and Sloan wgrriz-bands in both the SDSS and 2SLAQ cat- 
alogues. Richards et al. (2005) compared the h,j and g pass 




4000 5000 6000 7000 

Wavelength (A observed frame) 



8000 



Figure 1. Example of a 2dF spectrum with poor response at 
the blue end. In the top panel the original spectrum is shown as 
reduced by 2dFDR, the dashed line shows the response correction 
from Lewis et al. (2002) scaled for comparison. In the bottom 
panel this correction is applied to the spectrum. We find that, 
while this does improve the shape of the spectrum, it is inadequate 
in this case to correct for the sharp drop in received flux at the 
blue end. 



bands and found them to be roughly equivalent. Examin- 
ing 2QZ QSOs with g-band imaging they found a consistent 
{g — bj) — —0.045. In our analysis we use p-band magnitudes 
for all objects taken from the SDSS and 2SLAQ datasets, 
and for the 2QZ we apply this correction to their bj magni- 
tudes. 

Throughout this paper we use a flat (nm,f^A) = 
(0.3,0.7), Ho = 70kms~^ Mpc~^ cosmology when calculat- 
ing absolute magnitudes. We use the K-corrections laid out 
by Gristiani & Vio (1990) and correct for galactic extinction 
as advised in the relevant catalogue paper. 



4.2 Line measurements 

Our spectral analysis routine implements two separate fit- 
ting processes which are linked via iteration. Firstly a com- 
bined iron and continuum model is fit to the region of the 
spectrum not affected by the Mg ll line, then a Gaussian pro- 
file is fit to the line itself. We iterate the procedure to more 
accurately define the region of the spectrum affected by Mg ll 
emission, and so improve our iron and continuum fits. This 
significantly improves our spectral line calculations. 
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4-2.1 Continuum and iron correction 

Continuum emission from QSOs is well described by a power 
law in the optical-UV region. However, 2dF spectra are not 
flux calibrated. Variations in response modulate the shape 
of 2dF spectra most often characterised by a drop oS in flux 
at the blue end. Since we cannot accurately model these 
spectra with a power law we describe the continuum with 
a quadratic fit. This can then simultaneously correct for re- 
sponse effects in 2dF spectra and, since we only flt this to 
a limited region of the spectrum, approximate a power-law 
shape for Sloan objects. Fig.2 shows our iron and continuum 
flts to both a 2dF spectrum which has a drop off in flux at 
the blue end, and a well-calibrated SDSS spectrum. 

To remove iron emission from our spectra we use the 
template for QSO iron emission derived from the narrow-line 
Seyfert I object I Zw 1 by Vestergaard & Wilkes (2001). This 
template is derived from an intrinsically narrow-line object 
which made it possible to isolate the species responsible for 
each component of the emission. Unfortunately, because the 
template is derived from a real object, it is least well defined 
in the region around strong emission lines where de-blending 
differing species becomes difficult. 

This limitation is of particular importance directly be- 
neath the Mg II line where the template shows no iron emis- 
sion. This is not due to a quantified lack of emission, but to 
the difficulty of de-blending the iron and Mg ll emission in 
this region. It has now become common practise to add flux 
to the template in this region (e.g. Kurk et al. 2007; Salvian- 
der et al. 2007) to make it more consistent with theoretical 
models (Sigut & Pradhan 2003). We follow the method of 
Kurk et al. (2007) and add a constant level of flux under 
the Mg II line at 20 % the intensity of the average in the 
2930 - 2970 A (rest frame) region. We flnd that the level 
of flux included under the Mgll line does affect our results. 
Given the range of likely levels of iron emission in this region, 
however, the affect is not significant and does not impact our 
conclusions. This will be discussed further in section 5.2. 

Since the iron template is derived from a narrow-line ob- 
ject it must be smoothed to properly describe the iron emis- 
sion in a broader- line QSO. We follow the same procedure as 
previous authors (e.g. Boroson & Green 1992) and make a 
selection of iron templates smoothed by Gaussians of width 
500, 750, 1000, 1500, 2500, 5000, 7500 and 10,000 km/s. We 
fit all of these templates to this spectrum and accept the 
best fit in terms of the x^ ■ 

We fit for continuum and iron emission simultaneously 
with the SVDFIT routine (Press et al. 1989). The fit is per- 
formed on a region bounded at the blue end at 2450 A by 
the [Ne iv] A2424 line, and at the red end at 3085 A which 
is the limit of the iron template (all wavelengths are rest 
frame). We also mask the Mgll line out of our fitting. 

On the first attempt we use a fiducial ±50 A mask for 
the line but on subsequent iterations we use a Gaussian fit 
to the line to define its boundaries. 

It is important to keep in mind that once we have sub- 
tracted the iron and continuum contributions from the data 
we introduce a strong covariance into our spectrum and in- 
dividual pixel values can no longer be thought of as inde- 
pendent in subsequent error calculations. 



4-2.2 Gaussian fitting, iteration and line width 
calculations 

We subtract the iron and continuum contribution to the 
spectrum and fit a single Gaussian to the Mg ll line with the 
MRQMIN routine (Press et al. 1989). We then perform the 
iron and continuum fit to the original spectrum again mask- 
ing out data within ±1.5FWHM of the centre of the fitted 
Gaussian. This process is repeated until successive Gaussian 
fits have FWHMs consistent to within 1/2 their error. 

In this analysis some thought must be given to what 
statistic to use when defining the width of a spectral line. 
The width of single and/or multiple Gaussian/Lorentzian 
fits to the line are often used to describe the line width 
(e.g. McLure & Dunlop 2004; Shen et al. 2007; Wang, Lu 
& Zhou 1998). We have already performed Gaussian fits to 
our line and we do find that despite the clear non- Gaussian 
nature of the Mgll line these give a reliable measure for 
the width when compared with the other methods outlined 
here. However, we find non-parametric estimators for the 
line width more appealing for this analysis, primarily be- 
cause they make no assumptions as to the line profile, but 
also because their errors can readily be calculated including 
the contribution of covariance introduced into the spectrum 
in our iron and continuum correction. 

A common statistic used for describing line width is the 
FWHM (Vestergaard 2002; Vestergaard & Peterson 2006). 
In high S/N spectra this gives a good determination of the 
line width, and the error for the FWHM can readily be cal- 
culated including the contribution from covariance in the 
spectrum. Unfortunately the measured FWHM is quite sus- 
ceptible to noise. How to define the maximum flux density of 
a line in an unbiased way is unclear and multiple crossings of 
the half maximum value in noisy spectra demand some sort 
of averaging which can also bias the measurement. Further 
we flnd that the susceptibility of the FWHM to the values 
of a relatively small number of pixels in a spectrum makes it 
an unstable and often inaccurate measure of the line's true 
width. 

Another measure of the line width that is becom- 
ing increasingly widespread is the dispersion of the line a 
(Fromerth & Melia 2000; Vestergaard & Peterson 2006; Wil- 
hite et al. 2007). However, we flnd that the excessive weight- 
ing this measure gives to the values of pixels in the wings 
of the lines makes it an unreliable estimator of line width in 
low S/N spectra. 

Finally inter-percentile values (IPV) can be measured 
(Whittle 1985). While at first glance the process of mea- 
suring an IPV width is similar to measuring the FWHM, 
the dependence of IPV widths on the cumulative flux distri- 
bution rather than the flux density at a given point makes 
the IPV measurements considerably more robust. Like the 
dispersion, IPV widths are somewhat affected by noise in 
the wings of lines; in particular this can affect the total flux 
of a line and how one deflnes the zero-point of the cumula- 
tive flux distribution. However, when calculating the disper- 
sion the weight given to a single pixel is proportional to the 
square of the displacement of that pixel from the line centre. 
This power-of-two dependence makes a very susceptible to 
noise in the wings of a line which is not a problem for IPV 
widths. 

Our calculation for the IPV width is performed only on 
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Figure 2. These figures illustrate the iron and continuum fitting process on two very different spectra. In each case the top panels show 
the Mgll line region in the original spectrum, heavy dashed and solid lines show the continuum and iron + continuum model fit to it 
respectively. The bottom panel then shows the residual when this is subtracted from the spectrum, (a) shows the fit to the spectrum of 
J095421. 6-000152 as observed as part of the 2QZ. This is the same spectrum as given in Fig.l and shows a strong turnover in received 
flux at the blue end. (6) shows the spectrum of J000042. 89+005539. 5 as observed with the Sloan telescope and exhibits the classic quasar 
power law continuum. For both cases the quadratic continuum we use produces an excellent fit to the data. In all plots the y axis gives 
the normalised flux. 



the part of the spectrum within ±1.^FWHM of the Gaus- 
sian fit to the Mg II line. In this region we find the cumulative 
flux distribution and search this for the first crossing points 
at 1/4 and 3/4 the total flux. We then interpolate between 
the adjacent pixels either side of the crossing to obtain sub- 
pixel accuracy in our IPV estimate. 

Errors on the IPV widths are calculated from the spec- 
tral variance array including the contribution of covariance 
introduced by the iron and continuum subtraction. We find 
the mean error of these measurements to be ±0.05 dex. 

Finally we subtract the resolution of the spectrograph 
in quadrature from our measured line width under the as- 
sumption of a Gaussian profile for both the emission line and 
instrumental resolution. Since the Mgll line is not Gaussian 
in shape this is only an approximate correction and to test its 
validity we evaluate the effect of smoothing Lorentzian pro- 
files by a Gaussian. In the worst case scenario of a Lorentzian 
line with the smallest IPV width we measure in our sample 
(~ 1200 km/s; see Fig. 6), smoothed by a Gaussian with 
FWHM 675 km/s (equivalent to the 9 A resolution of the 
2dF spectrograph at A = 4000 A) the correction above re- 
sults in an overestimate of the line width by < 5 %. For 
broader spectral lines, and in cases where the spectral reso- 
lution is greater, this effect is diminished and hence we make 
no further attempt to correct for instrumental smoothing. 
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Figure 3. Comparison of 50% IPV widths anfd FWHMs for 
our sample. The dashed line gives the 1:1 relation and the solid 
line shows the y on x least squares flt. IPV widths have been 
multiplied by 1.75 for the comparison which is the expected ratio 
for a Gaussian line 



4.2.3 IPV widths vs. F WHM 

In taking the IPV width as our measure of line width we 
assume there is a linear scaling relation between the differing 
line width parameterisations. Specifically we assume that 
equation 3 holds equally true for IPV widths. To test this 
assumption we also measure the FWHM of the Mg ll lines 
in our sample to compare with the IPV widths. 

These are measured by defining the maximum flux den- 
sity of the spectral line as the average of the highest pixel 



value and the two pixels adjacent to it. Then the lower and 
upper bounds of the FWHM are found by averaging all 
crossings of the half maximum value on the red and blue 
sides of the lines center. Fig. 3 shows a comparison between 
the IPV width and FWHM measurements for our dataset. 
In the figure we have multiplied the IPV widths by 1.75, 
the ratio between IPV and FWHM for a Gaussian. There 
is considerably more scatter in the FWHM direction due to 
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the susceptibility of this statistic to noise in the spectrum. 
However, a y on x least squares fit to this data gives 



\og{FWHM) = 0.981og(JPV) + 0.1, 



(4) 



a linear relation between the IPV width and FWHM. The 
offset is due due to the profile of the Mg ll line which tends 
to have more flux in the wings of the line compared to a 
pure Gaussian. 



4.2.4 Rejection eriteria 

With an automated routine it is impossible to accurately 
fit every spectrum and there is always the chance of catas- 
trophic failure. We make a series of checks and cuts to the 
results to try and reduce the number of spurious fits in our 
final sample. 

Failure to produce a satisfactory fit is normally due to 
a property of the spectrum being analysed. Absorption fea- 
tures, low S/N, extreme curvature in the underlying spec- 
trum, bad pixels and residual telluric features can all cause 
failures, and commonly a failure is due to a combination of 
these. 

To test the effect of S/N on the accuracy of our fit- 
ting routine we simulate large numbers of noisy spectra by 
adding Gaussian noise to high quality Sloan and 2dF spec- 
tra. We then measured the width of the Mg ll line in these 
degraded spectra and compared our results with the mea- 
surement from the high S/N spectrum. 

Down to a S/N ~ 1.5 pix^^ wc find no systematic devi- 
ation between the average line width measured from noisy 
spectra and the true line width. Furthermore we find that 
above this S/N the errors on our measurements correctly 
describe the scatter we observe in the measured line widths. 
We therefore apply a spectral S/N cut of S/N > 1.5 pix^^ 
to our sample. 

It is worth noting that while our error calculations as- 
sume Gaussianity in linear space we find that the distribu- 
tion of IPV widths measured for a given spectrum at low 
S/N is better described by a log-normal. This may not be 
surprising since the widths are randomly distributed and 
limited to be > 0, however, that our errors transfer from 
linear to log space is encouraging. We transform the linear 
error aun to the log error oiog via its ratio to the measured 
line width IPV . If we propagate the error through the log 
then 



2 _ 
f^log — 



d 



Clog 



dIPV 

O'lin 



(log,„(JPV))) aL 



/py In(lO)' 



(5) 
(6) 



We make three further constraints on our data as to 
whether we will accept a particular fit. We apply a redshift 
cut to avoid contamination by the many telluric features to- 
wards the red end of our spectra, we limit the number of 
times we iterate our fitting procedure as described in sec- 
tion 4.2.2, and finally we try to eliminate broad-absorption 
line (BAL) objects which could contaminate our data. 

Telluric features pervade unreduced spectra redwards of 
7100 A. While both SPECTR02d and 2dFDR try to remove 
sky features the reduced spectra often contain residuals in 
this region. To avoid the worst of this we make a simple 




Figure 4. Distribution of our final sample of objects in the 
redshift-luminosity plane. Magnitude limits for the various sur- 
veys can be made out as lines where the density of objects in- 
creases. The very few objects with g-band magnitudes below the 
21.85 cut for 2SLAQ were selected as potential high redshift ob- 
jects in the i-band. 



redshift cut at 2 = 1.5 which ensures that the Mgll line and 
much of our continuum fit avoid this region. 

For a typical spectrum we iterate the fitting procedure 
outlined above 2-3 times to obtain convergence. We limit 
the maximum number of iterations allowed to 20 and reject 
any fits which have not converged by this time as unreliable. 

BAL objects are a contaminant in our data and we try 
to identify and reject these during our fitting. In the 2QZ 
and 2SLAQ catalogues many of the most severe BAL objects 
are flagged as such and are not included in our analysis. 
Those which are not flagged in a catalogue are identified on 
the basis that they will have pixels which are significantly 
deviant from the rest of the spectrum. After every Gaussian 
fit we reject any pixels which deviate by > 3a from the fit 
and we do not include these pixels in subsequent iterations. 
Once convergence is obtained in our fitting process we reject 
any object which has consecutive rejected pixels spanning 
500kms"\ 

Note that this process of rejecting pixels does not only 
affect BAL objects; any pixels with outlying values will be 
excised. If this process removes more than one fifth of the 
pixels within ±1.5FWHM of the Mgll line we also reject 
the fit from our data. 

Details of how many objects are rejected due to these 
various cuts are given in table 2. In total ~ 4 % of the objects 
are rejected from our sample. More than 80 % of these are 
objects rejected due to the fitting routine not converging on 
a solution which occurs primarily on low S/N spectra. This, 
combined with the imposed S/N cut could create a bias at 
low luminosity in our data, in particular if there is a strong 
correlation between IPV width and luminosity. However, as 
we shall see there is only slight evidence for a correlation 
between luminosity ind IPV width. In addition to this the 
small number of objects rejected mean that any bias will be 
negligible. On the other hand the BAL rejection is designed 
to expel true outliers from our data and should not bias the 
results in any way. 
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Table 2. Table outlining our final dataset. For each of the samples wc give the total number of objects for which we have attempted to 
analyse the Mg II line as well as the apparent (uncorrected for Galactic extinction) and absolute (extinction corrected) magnitude range 
of these objects. We also give the number of these fits which have been rejected from our analysis by each of our criteria. These are low 
S/N spectra, possible BAL objects, and objects for which the fitting code did not converge on a result within twenty iterations. Note 
that the redshift cut at z = 1.5 to avoid sky contamination has already been factored into the number of spectra in column 2. We also 
give the total number of fits rejected, since an object can be rejected for a number of reasons this does not equal the sum of the previous 
3 columns. 













No. of fits rejected 




Sample 


No. fit 


App. Mag. Range 


Abs. Mag. range 


S/N 


BAL 


ltt.>20 


All 


2SLAQ 


2684 


25.09 > g> 18.12 


-19.48 > Mg > -26.38 


44 


53 


142 


193 


2QZ 


7209 


20.80 > g> 18.20 


-21.07 > Mg> -26.74 


161 


75 


356 


425 


SDSS 


23725 


23.33 > g> 15.32 


-20.98 > Mg> -29.27 


83 


141 


653 


786 


All 


33618 


25.09 > g> 15.32 


-19.48 > Mg> -29.27 


288 


269 


1151 


1404 




Figure 5. This plot gives the results of our line fitting analysis. We plot the absolute g-band magnitude of the source vs. the measured 
IPV velocity width of the Mgll line. Over plotted are lines of constant SMBH mass (dotted) and Eddington ratio (L/I/Bdd; dashed) as 
a guide to where these objects fall in mass-accretion space. Masses are labeled in units of Mq. We have added contours of even density 
to the plot to highlight the shape of the distribution. The contours are evenly spaced in terms of log(density). 



Fig. 5 shows the results from our fitting plotted on a line 
width - magnitude diagram, and Fig. 4 shows how these ob- 
jects are distributed on the redshift- luminosity plane. Each 
point represents an object with an accepted line width mea- 
surement. The fitting results for these objects is available 
from the 2SLAQ website (www.2slaq.info). 

Added to the plot are lines of constant SMBH mass 
(dotted) and Eddington ratio (dashed). These are calculated 



using the McLure & Dunlop (2004) calibration for the Mg ll 
virial relation, assuming their B-band bolometric correction 
and taking g ~ B — 0.11 to calculate the bolometric lu- 
minosity of the sources. We then use a bolometric correc- 
tion to the continuum luminosity at 3000 A of 5.2 (Richards 
et al. 2006) to calculate the continuum luminosity at this 
wavelength from the bolometric luminosity. The McLure & 
Dunlop calibration was based on line widths measured in a 
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1000 2000 5000 

IPV Velocity Width (km/s) 

Figure 6. Distribution of Mgll line widths for our sample (heavy 
histogram) with a Gaussian fit for comparison (fine line). We find 
this to be roughly log-normal although with slightly more objects 
in the wings in particular towards larger line widths. 



different fashion to this work. They fit two Gaussians to the 
Mg II fine and take the FWHM of the broader component as 
their hne width. We correct the lines in Fig. 5 for our use of 
IPV widths assuming a Gaussian profile for the Mgll lines, 
but McLure & Dunlop's use of two Gaussians in their fitting 
will likely bias the plotted lines high by a factor of 1.5 to 2. 
These lines are plotted more as a guide as to how lines of 
constant mass and Eddington ratio would lie on the diagram 
rather than being exact in normalisation. 

We will discuss the implications of this diagram further 
in section 6, but already if one compares contours on the left 
and right hand side of this plot it appears to suggest that 
there is less scatter in broad-line widths for more luminous 
quasars. Furthermore the average line width shows very lit- 
tle variation across the luminosity range sampled, a result 
consistent with previous analyses (Corbett et al. 2003; Shen 
et al. 2007). 



4.3 Dispersion analysis 

To investigate how the dispersion in Mbh depends on the 
luminosity of QSOs in our sample we bin our sample by 
luminosity and calculate the dispersion in line widths in each 
bin. To remove possible redshift evolution in this relation 
we also bin by redshift. As with measuring line widths some 
consideration must be given to choosing the most robust 
estimate for the dispersion in each bin. 

The overall distribution of QSO line widths is shown 
in Fig. 6. We find it is roughly log-normal although with ex- 
tended wings and a slight asymmetry. This further motivates 
our use of logarithms in equation 3 to derive the dispersion 
in Mbh. 

The obvious way to measure dispersion is with the rms 
in log space. This approach is, however, very susceptible to 
outliers. Objects with either very large or very small results 
for their line width are more likely to be due to poor spectral 
fitting and tend to have large errors on these measurements. 
We find that the IPV method is less biased by these out- 
liers although for the most part we find very little difference 
between the two statistics. 

We must be wary that the distribution of measured line 



widths is not a direct representation of the intrinsic distribu- 
tion of line widths in our sample. Instead the measured dis- 
tribution represents the intrinsic distribution convolved with 
the distribution of errors on these estimates. In the perfect 
case where both are normally distributed in log space the 
intrinsic rms, rmsi, is related to the measured rms, rmsm, 
and the average error in line width a^ by 



= vr 



(7) 



We are not using the rms in this analysis, instead we 
take the 68.3 % IPV width to provide an equivalent statistic. 
Further we do not take a^ to describe the average error in 
our data since we find this can be skewed by a small number 
of objects with very large errors on their line widths. To 
avoid these outliers we take the square of the median error 
in a bin to estimate ct^. 



5 RESULTS 

We calculate the dispersion in line widths for each L — z bin 
and correct by the average error on the line widths in that 
bin. To ensure reliable dispersions we only consider L — z 
bins that contain more than 40 objects. Fig. 7 displays this 
corrected dispersion vs luminosity for each of the 2QZ (a), 
2SLAQ (6) and SDSS (c) samples separately, and then for 
the combined dataset (d). Note that while we bin our data 
by redshift we do not see evidence for redshift dependence 
in our sample (see section 5.1 below) and so we plot the data 
from all redshift bins in Fig. 7. 

Error bars in the plots are propagated from the errors 
on the line width measurements. For the most part we have 
a large number of objects in each L ~ z bin and thus are 
dominated not by these random errors but by systematics 
in the analysis. 

Figs. 7(a) and (6) have considerable scatter and exhibit 
little-to-no trend in the plots, although the shallow trend 
seen in the 2QZ data is significant (a Spearman rank test 
gives rs — 0.56 with P{rs) = 0.01). Fig. 7(c) shows a strong 
relation with little scatter due to the larger numbers in the 
SDSS sample (Table 2) and its higher quality spectra. 

Fig.7(rf) shows the results when we combine all of our 
data. This is heavily dominated by SDSS objects down to 
absolute magnitudes Mg ~ — 23 where the 2QZ and then 
2SLAQ samples become important. 

Shen et al. (2007) perform a similar analysis on the 
SDSS DR5 quasar sample. They find a consistent level of 
variance in their data but report no dependence on luminos- 
ity, although inspection of their Fig. 4 does suggest a trend 
albeit slight. It may be that their following the McLure & 
Dunlop (2004) prescription for measuring the width of the 
Mg II line and/or the use of large luminosity bins in this fig- 
ure could explain their not finding as strong a relation as 
apparent in Fig. 7. 



5.1 Redshift dependence 

We find in Fig. 7 that points from the same luminosity bin 
but differing redshift bins lie on top of each other in the 
plot suggesting the dispersion in line widths has little-to-no 
redshift dependence. This is further illustrated in Fig. 8. In 
this plot we show the dispersion in IPV line width plotted 
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Figure 7. Plots of the dispersion in log(/PV') in absolute magnitude bins for the 2SLAQ (a), 2QZ (6) and SDSS (c) samples, as well as 
the combined sample (d). Since we find no redshift dependence in our data we plot the results from all redshift bins on the one figure, 
hence more than one point per magnitude interval. The open symbols show the dispersion in measured line width before correction 
{rmSm', see equation 7) for errors on our measurement, filled symbols with error bars give the corrected dispersion (rmsi). 



against redshift for L — z bins. Note that these are not the 
same L — z bins as used in the rest of this paper, we have 
doubled the number of redshift bins and halved the number 
of luminosity bins to make the plot clearer. 

In the plot we connect points which are in the same 
luminosity bin and it is evident the dispersion in IPV line 
widths in these bins has at most a weak redshift dependence. 
A Spearman rank test performed only on the three luminos- 
ity bins which have full redshift coverage gives rs = —0.14 
significant at only the 50 % level. This is consistent with the 
changing distribution of luminosities within each L bin as 
the luminosity function evolves with redshift. 



5.2 The effects of the iron template 

As discussed above, the Vestergaard & Wilkes (2001) iron 
template used in our analysis includes no iron emission di- 
rectly under the Mg ll line and we correct for this by adding 
flux at a level suggested by theoretical models for QSO iron 
emission (Sigut & Pradhan 2003). To test what effect this 
has on our results we refit all of our data with an iron tem- 
plate with no iron emission under Mgll and twice as much 
iron emission as used in the main analysis. Finally we also 
try an iron template where we interpolate between the iron 
emission peaks at ~ 2750 and 2840 A (see e.g. Fig. 13 Sigut 




Figure 8. Here we plot the dispersion in log(7Py) vs. redshift 
for our sample. Luminosity bins are 1 mag. in width and are cen- 



tered on (top to bottom) Mg 



-21, -22, -23, -24, -25, 



-26 



and —27. Bins with the same luminosity cuts are given the same 
symbol and connected in the plot, and it is evident the dispersion 
in IPV widths in these bins is only very weakly dependent on 
redshift. 
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& Pradhan 2003); this template includes more than five 
times the flux in the Mgll region as in our primary tem- 
plate. 

There is a systematic trend between our results and the 
amount of iron emission added under the Mgll line. The 
more emission, the more scatter we find in our line widths. 
However, the effect is small. We find an offset of ~ 0.015 
between the dispersion in log(/PV') for objects fit with no 
iron under the Mgll line and twice as much iron emission 
as assumed in this analysis, the data presented in Fig. 7 
lie in between these two. In the case of the template where 
we interpolate between the iron peaks either side of Mgll 
we find an offset of ~ 0.08. While this approach produces a 
significant offset this model has no physical basis and was 
only implemented to test the extremes of the distribution. 

Furthermore, regardless of which iron template we use 
in our fitting the trend in our dispersion analysis is the same. 
We always find that more luminous objects show less scatter 
in line widths than fainter objects. 
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Figure 9. This shows how the skew in the IPV line width dis- 
tribution is dependent on the magnitude of QSOs. The skew is 
defined by equation 8 and the error bars are calculated from the 
errors on the IPV width measurements propagated through this 
equation. 



5.3 Completeness and homogeneity 

While the 2SLAQ, 2QZ and primary SDSS QSO surveys 
all have high spectroscopic completeness overall (~ 90%), 
towards the faint flux limit of each survey this drops off. 
Furthermore, the amalgam of many separate surveys in the 
SDSS quasar catalogue and our adding to this the 2SLAQ 
and 2QZ spectra as well makes for a very inhomogeneous 
sample of objects in terms of selection criteria. However, 
these factors are not of great importance in terms of biasing 
our results or the trend shown in Fig. 7. 

Incompleteness at faint magnitudes should not have a 
major affect since inspection of Fig. 5 shows there is little 
variation in the average line width with luminosity. Hence 
the dispersion in IPV widths in a luminosity bin is not af- 
fected strongly by the size of a luminosity bin, or the distri- 
bution of objects within a given bin. 

Inhomogeneity in selection criteria is also not respon- 
sible for the observed trend in Fig. 7. This is evident since 
we find exactly the same trend in each sample separately. 
Furthermore we find the same trend if we only take the pri- 
mary SDSS quasar sample, albeit with significantly reduced 
luminosity coverage. 



6 DISCUSSION 

In Fig. 5 the results of our line fitting are plotted along 
with contours of constant SMBH mass and Eddington ra- 
tio. While there is uncertainty in the normalisation of these 
lines the distribution of points in the plot relative to them 
is quite suggestive. The line width distribution appears to 
be constrained on at least two sides parallel to these lines. 

Firstly the bottom of the distribution follows the line at 
L/L-Edd ~ 1- Potentially this indicates that the Eddington 
rate does represent an upper limit for allowable accretion in 
QSOs, constraining our results to lie above this line. 

Secondly brighter than Mg ~ —24.5 the top of the dis- 
tribution of line widths follows the Mbh = 10^° M© contour 
on the plot. This is likely due to the drop off in the SMBH 
mass function at high mass. The space density of SMBHs 



falls off dramatically for masses above 10® Mq and this cre- 
ates an upper limit to the distribution of line widths plotted 
in Fig. 5. 

Finally, and more speculatively, fainter than Mg ~ 
—24.5 the top of the line width distribution appears to be 
constrained along a line of constant Eddington ratio around 
L/L-Edd ~ 0.01. The implication of this drop off of objects 
towards low accretion rates is that there is a preferred level 
of accretion onto SMBHs for QSOs. This indicates that the 
accretion distribution for SMBHs as a whole is bimodal, 
with quiescent and active SMBHs occupying distinct areas 
of mass-accretion space. 

As a test as to how the Eddington limit effects the ac- 
cretion efficiency distribution we measure the skew of the 
line width distribution in magnitude bins. We take the di- 
mensionless skew defined by the third moment of the IPV 
width distribution 



skew ■ 



1 V^ / IPV^ - IPV 



N ^-^ 

i = l 



aipv 



(8) 



Fig. 9 plots the skew in the IPV width distributions against 
the magnitude of the bins. 

This figure shows that in more luminous QSOs the line 
width distribution consistently has a positive skew. This 
skew towards larger line widths and hence lower accretion 
efficiencies may be the result of the accretion efficiency dis- 
tribution being truncated at the fast end by the Eddington 
rate. 

There is some evidence that this skew reverses for 
fainter objects. At these luminosities the mean of the ac- 
cretion efficiency distribution is significantly lower than for 
brighter objects, hence a cut at the Eddington rate will not 
have as strong an effect. The reversal of the trend is likely 
due to the underlying mass function of active SMBHs. With 
many more small SMBHs there would be a skew towards 
higher accretion efficiencies and hence smaller broad line 
widths at a constant luminosity. 

Taken at face value the distribution of points in Fig. 5 
fits well with simple expectations about the QSO SMBH 
population, suggesting that the virial mass estimates work 
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relatively well. However, see section 6.3.1 for a discussion on 
some concerns with the virial mass technique which these 
data highlight. None the less, this gives us confidence that 
we can make a meaningful comparison between our data and 
theoretical models for activity in the SMBH population. 



6.1 Comparisons with models 

In Fig. 10 we compare our results to the model predictions 
of Hopkins et al. (2005b). The models are dependent on 
the distributions of galaxy merger rates at a given epoch. 
This function is not well constrained by current observations 
and two possible realisations are applied. Parametrically the 
merger rate function is described by a double power law and 
the differing models assume differing slopes for the low mass 
end. The first assumes a steep drop off of merging systems 
towards lower mass and gives only a narrow mass range for 
merging systems at a given redshift. Secondly a broad range 
of systems are assumed to be merging at any one time with 
a fiat slope at the low mass end of the merger function. We 
label these models as sharp or flat as describes their low 
mass slope; the sharp and fiat models are shown as the solid 
and dashed lines respectively in Fig. 10. 

We find that our data points match the Hopkins et 
al. (2005b) models, extremely well in general. However, the 
models exhibit some evolution with redshift which we do not 
observe in our data. 

These models are, to an extent, bound to the QSO lumi- 
nosity function which shows strong redshift evolution. Over 
the redshift range we sample the break in the luminosity 
function, M* , changes by more than a magnitude, and the 
Hopkins et al. models follow this to an extent. We find no 
such trend in the relation shown in Figs. 7 and 10. To illus- 
trate this Fig. 11 shows exactly the same data as Figs. 7(d) 
except the dispersion is plotted against M* — M, i.e. the 
object's luminosity relative to the break in the luminosity 
function. For this plot we take the quadratic parametrisa- 
tion of M* from Groom et al. (2004) . Comparisons between 
Figs.7((i) and 11 clearly show that the observed trend to- 
wards lower dispersion in broad-line width for brighter ob- 
jects is dependent on the object's luminosity, not its position 
on the luminosity function. 

The Hopkins et al. models evolve with redshift in a simi- 
lar sense to the QSO luminosity function albeit less strongly, 
hence we observe a slight discrepancy. However, it was un- 
likely that these simple models would be exact in their de- 
termination of the dispersion in Mbh. Of more importance is 
the trend, which is predicted by the models and well echoed 
in the data. Our data show that more luminous QSOs show 
less scatter in their Mg ll line widths than those of lower lu- 
minosity, implying they have less scatter in their black hole 
masses. 



6.2 Is this trend a selection effect? 

Recently Babic et al. (2007) showed that a broad intrin- 
sic distribution of Eddington ratios which truncates at the 
Eddington rate, convolved with a double power law SMBH 
mass function for QSOs, naturally leads to a selection effect 
such that samples with a fainter flux limit will find a broader 
range of Eddington ratios. 
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Figure 11. Same data as Fig. 7(d), except the dispersion is plot- 
ted vs. M* — M rather than M. Comparing the scatter evident in 
this plot with Fig. 7(d) clearly demonstrates the observed trend is 
dependent on the brightness of a quasar, not its position on the 
luminosity function. In this plot we use the quadratic parametri- 
sation for M* given in Croorn et al. (2004). 



To test whether this effect could be responsible for the 
trend observed in Fig. 7 we recreate this situation and com- 
pare it with our results. In our analysis we assume that the 
distribution of Eddington rates (L/I/Edd ~ A) is log-normal 
with a given mean, (A), and dispersion, ax, and is truncated 
at A = 1. To begin with we assume the SMBH mass function 
can be described as a double power law, and constrain this 
by convolving with the Eddington distribution and fitting 
to the observed luminosity function (we use the Hopkins, 
Richards & Hernquist 2007 B-band luminosity function at 
redshift 1), producing an (almost) unique solution for the 
mass function. We can then derive the dispersion in SMBH 
mass at a given luminosity. 

We find that the effect discussed by Babic et al. (2007) 
is reproduced in the situation where (A) and ax conspire 
such that the A = 1 cut is not too many ax 's away from (A) ; 
i.e. in situations where the cut has a significant effect on 
the Eddington distribution. Hence in situations where (A) is 
small ax must be large to produce the effect. And if ax is 
small then (A) must be comparable to 1 for a pronounced 
effect. Fig 12 gives an example of this. We plot the dispersion 
in SMBH mass in these models as a function of luminosity 
in the case where (A) = 0.3 and 0.7, for a variety of values 
of CTa. 

While both of these plots show a decrease in the dis- 
persion in SMBH mass with luminosity, neither give a good 
representation of our results. In the case where (A) = 0.3 ax 
must be so high to get a pronounced gradient that the nor- 
malisation of the relation, in particular at the bright end, is 
too high when compared with our results. In the case where 
(A) = 0.7, ax can be smaller and the normalisation is closer 
to that observed. However, the model dispersions fiatten a 
magnitude brighter than M* , an effect we do not observe. 

In all cases the dispersion asymptotes to a constant 
value above M* since at these luminosities one is sampling 
from both a luminosity function, and by extension SMBH 
mass function, which follow a single power law. The Babic 
et al. effect does not occur without sampling objects from a 
mass function which is steepening in the log sense. 
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Figure 10. Comparison between our data and the models of Hopkins et al. (2005b). The model lines give the predicted dispersion 
in log(MBH) as a function of luminosity at the mean redshift of the bin (top right of each panel). The two different curves give the 
predictions from differing merger rate functions, either with a sharp cutoff at lower masses (solid) or a flatter cutoff (dashed). Our 
measured dispersion in log(/PV') is multiplied by 2 for the comparison (equation 3). 
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Figure 12. These plots illustrate the effect described by Babic et al. (2007). The dispersion in SMBH mass at fixed luminosity is plotted 
assuming Eddington ratio distributions with means (A) = 0.3 and 0.7, and with a variety of widths (ux). 



This is, however, an artifact of our assigning a dou- 
ble power law for SMBH mass function. If we assume a 
Schechter function instead the high mass end drops off expo- 
nentially and the Babic et al. effect is apparent at the higher 
luminosities sampled. On the other hand with this model the 



fit to the observed QSO luminosity function (which does fol- 
low a double powerlaw) is significantly degraded. 

Finally, since these models are tied directly to the lumi- 
nosity function (i.e. are relative to M*) they evolve strongly 
with redshift. The observed evolution in the QSO luminos- 
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ity function must be either due to a similarly evolving active 
SMBH mass function or an equivalently varying accretion ef- 
ficiency distribution. Either will strongly affect the Babic et 
al. effect. 

As discussed in section 5.1 we find no evidence for a red- 
shift dependence in our results. And as shown in Fig. 11 the 
observed correlation between dispersion in IPV widths and 
luminosity is degraded when plotted vs. luminosity relative 
toM*. 

In light of the above discussion it is clear that while 
the Babic et al. (2007) effect may well bias our results, it 
does not accurately describe both the results presented in 
Fig. 7 and the observed luminosity function simultaneously. 
The lack of observed evolution in our results may suggest 
that this effect is not significant, however, without a better 
understanding of the underlying SMBH mass function or 
accretion efficiency distribution it is difficult to gauge the 
magnitude of this effect and we cannot rule out a significant 
bias in our data. 



6.3 Disp(log(/PV)) vs Disp(log(MBH)) and virial 
masses revisited 

While we observe a strikingly good agreement between our 
data and the the Hopkins et al. (2005b) we caution against 
reading too much into the normalisation of our data with 
respect to the dispersion in SMBH mass. In fact a closer 
examination of this normalisation raises a question as to 
just how accurately the dispersion in \og{IPV) traces that 
in log(MBH). 

In using equation 3 we are assuming there are no other 
sources of scatter in the QSO SMBH population which will 
affect these results. This is almost certainly not the case and 
we must consider where other sources of scatter might oc- 
cur, their magnitude and their affect on our results. Shen et 
al. (2007) recently performed simulations of the QSO SMBH 
population to investigate the biases associated with virial 
mass estimations and we follow a similar path to their work. 

They describe the QSO and SMBH populations with 
four variables. The SMBH mass, the bolometric luminosity 
(and by extension the Eddington ratio), the monochromatic 
luminosity as used in the virial relation (equation 2) and 
the FWHM of the spectral line being analysed (in our case 
the IPV width of Mgll). We add to this the radius and 
velocity dispersion of the BLR: tblr and ublr. There are 
possible sources of scatter to each of these variables and here 
we discuss each of them, their likely magnitude, and their 
affect on our measurements. 

• Scatter in luminosity 

The primary source of scatter in luminosity at fixed 
SMBH mass comes from the distribution of accretion effi- 
ciencies for QSOs. This is precisely the scatter we are trying 
to measure in this work and is not a contaminant to our 
calculations. 

The only other sources of scatter to the luminosity are 
photometric errors, host galaxy light and extinction. Each of 
these would add extra scatter into the luminosities we mea- 
sure independently from any virial mass equations. Hence 
this scatter would not propagate through to the line widths 
we measure in the way scatter in accretion rate would. This 
has the unique effect of making the SMBH mass distribution 



we measure at a given luminosity narrower relative to the 
intrinsic mass distribution (see Shen et al. 2007). 

However, these sources of scatter can only play a mi- 
nor role in our measurements and calculations. Photomet- 
ric errors are ~ 0.1 mag. in the fej-band used by the 2QZ 
(~ 0.04 dex in L; ^ 0.02 dex in Mbh) and smaller for the 
SDSS observations. Host galaxy emission and extinction will 
only affect the fainter objects in our sample, and even at 
these luminosities will not be a major factor. Hence these 
sources of scatter will not be significant to our calculations. 

Shen et al. also include a scatter due to variations in con- 
tinuum shape when converting bolometric luminosities to 
monochromatic luminosities. They quote a value for this of 
~ 0.1 dex translating to ~ 0.05 dex scatter in Mbh and so 
again this is not of great importance to our analysis. 

• Scatter in rBLR 

The r — L relation shows that the radius of the BLR 
around QSOs is driven by the luminosity of the source. 
Hence any scatter in the intrinsic luminosity of the source 
(i.e. not due to photometric errors/host galaxy emis- 
sion/extinction) will also scatter tblr. Furthermore the r—L 
relation has been observed to exhibit intrinsic scatter of 
its own at a level of 40% (Kaspi et al. 2005) translating 
to at least 0.15 dex, and this intrinsic scatter will further 
broaden the tblr distribution. This extra scatter will prop- 
agate through to the IPV width distribution we measure, 
and will bias our final measurement for the scatter in Mbh 
high by the same factor. 0.15 dex of scatter in Mbh is signif- 
icant to our discussion and we expand on the implications 
of this in section 6.3.1. 

• Scatter in ublr 

The distribution of velocity dispersions will have all of 
the above scatters folded into it apart from the scatter 
in luminosity due to photometric errors/host galaxy emis- 
sion/extinction. If we take wblr to be the virial velocity of 
the BLR then there will be no additional sources of scatter 
which affect this quantity since the virial equation (eqn. 1) 
is exact. 

However, other (non-virial) factors may effect the velocity 
of the BLR. Shen et al. include an additional scatter to the 
line width distribution due to non-virial velocities in the 
BLR. This will bias the scatter we measure in IPV widths 
(and hence Mbh) high, although in this work we do not 
attempt to investigate this further. 

• Scatter in IPV widths 

Beyond all of the scatters which affect ublr, the IPV 
width distribution we measure can be effected by further 
sources of scatter. Firstly errors on our measurements will 
artificially broaden the measured IPV width distribution, 
and we account for this in our dispersion analysis (see sec- 
tion 4.3 and equation 7). 

Secondly the line width we measure may not accurately 
describe the velocity of the BLR. Selective absorption (e.g. 
the model proposed for the C iv BLR by Richards et al. 
(2002) and/or orientation dependent line widths would fur- 
ther broaden the measured IPV width distribution. There is 
little evidence for absorption in the Mgll line for non-BAL 
QSOs and do not attempt to model its effects here. The 
potential effects of orientation on our measurements are dis- 
cussed in section 6.4 below. 

Of these sources of scatter, only the uncorrelated scat- 
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ter in luminosity can bias the measured dispersion in IPV 
widths low. All of the other sources of scatter propagate 
through to the IPV distribution and will bias the results 
high. 

Little is known about the magnitude of many of the 
sources of scatter discussed above, but most can be ruled 
out as negligible when compared to the observed dispersion 
in IPV widths. The only source of scatter which is known to 
be comparable to that measured for the IPV widths is the 
intrinsic scatter in the r — L relation. 



6. 3. 1 Implications of intrinsic scatter in the r — L relation 

Equation 3 assumes that there is no variation in raLR within 
a luminosity bin. This is not the case, and the distribution of 
line widths in a given luminosity-redshift bin is not directly 
analogous to the intrinsic distribution of SMBH masses. In- 
stead the distribution of line widths is a convolution of the 
black hole mass distribution and the distribution of BLR 
radii in a given bin (as well as our error distribution; equa- 
tion 7). 

With this in mind a better estimate of the true disper- 
sion in Mbh is the calculated dispersion in \og{IPV) (x2; 
equation 3) with both the dispersion in luminosities in a 
single bin (xa; equation 2) and the intrinsic scatter in the 
r — L relation subtracted from it in quadrature. 

The scatter in log(L) for 0.5 mag. bins will be ~ 0.1 dex 
in the faint bins, and decrease to brighter luminosities as the 
luminosity function steepens. This translates to < 0.05 dex 
of scatter in Mbh and hence is relatively unimportant to 
this discussion. 

On the other hand the intrinsic scatter in the r — L 
relation is not negligible with respect to our results. Since 
the radius of the BLR is dependent on the luminosity of the 
quasar (and not the other way around), any intrinsic scatter 
in this relation will also tend to bias our results high. Due to 
the relatively small number of objects with good reverbera- 
tion data the scatter in the r — L relation is not definitively 
constrained, and not defined at all for the brighter luminosi- 
ties in our sample. However, Kaspi et al. (2005) do measure 
this for the H/3 BLR with 35 of the best mapped AGN and 
quote a value of ~ 40 % for the intrinsic scatter in tblr. This 
corresponds to at least ~ 0.15 dex in the log (0.17 if the er- 
ror is propagated as in equation 6), and is greater than the 
scatter we find in SMBH mass in the brightest magnitude 
bins. 

This leaves a dilemma: Either there is significantly less 
intrinsic scatter in the r — L relation than that quoted by 
Kaspi et al. (at least at high L), or there is an intrinsic 
problem with our approach to finding the scatter in Mbh. 

The 40 % scatter quoted by Kaspi et al. (2005) is not 
directly applicable to our analysis for two reasons: 1) The 
r — L relation studied in that work is for the H/3 emission 
region not Mgll and 2) The r — L relation is only defined 
for relatively faint Seyfert 1-type objects, while the lowest 
scatter we find is for the most luminous quasars in our sam- 
ple. In effect our use of the Kaspi et al. value for the scatter 
in this relation equates to an extrapolation over ~ 2 orders 
of magnitude in luminosity. 

In their analysis Kaspi et al. suggest possible sources 
of scatter in the r — L relation such as intrinsic reddening, 
contributions by the host galaxy, or effects of variability due 



to non contemporaneous observations. Of these, galaxy con- 
tamination and potentially reddening will play a smaller role 
at higher luminosities. Bentz et al. (2006) used HST imag- 
ing to correct a subset of the objects in Kaspi et al.'s sample 
for host galaxy emission and found that the intrinsic scatter 
in the r — L relation could be as low as ~ 30 %. At higher 
luminosities where host contamination would be lessened it 
is possible that the scatter is even less. 

If there were no scatter in the r — L relation our data 
imply there is only -^ 0.14 dex scatter in Mbh in luminous 
QSOs. This may be plausible considering that at the ex- 
treme luminosities we are considering here we would expect 
most quasars to be radiating at or around their Eddington 
luminosity. 

An alternative interpretation may be that the disper- 
sion in log(/PV^) does not properly represent the dispersion 
in log(MBH). That is, QSO broad lines do not show enough 
variation to account for the expected variation in Mbh when 
accounting for the intrinsic scatter in the r — L relation. 

This raises a question as to the validity of the analysis 
performed here, and the virial method for estimating Mbh. 

There is evidence for a virialised BLR. Time lags be- 
tween continuum and emission line variations have been ob- 
served to be faster in the wings as opposed to the core of 
lines (KoUatschny & Dietrich 1996; Onken & Peterson 2002; 
KoUatschny 2003). And the same authors also find that in 
objects which have had more than one line mapped their 
time delays appear to follow a virial (r oc V^^j^) relation. In 
addition, studies comparing virial masses to bulge velocity 
dispersion have shown a correlation analogous to that found 
in local quiescent galaxies (Onken et al. 2004; Woo et al. 
2006) indicating that the virial method works as a tracer of 
Mbh. 

An explanation could be that the BLR of QSOs is not 
wholly virialised. And there is a significant component of 
the BLR which shows very little object-to-object velocity 
variation. If this were the case virial motion would still be 
the primary cause of the variation between broad line widths 
in QSOs and observations such as the M — a relation would 
be reproduced. However, globally we would not see a range 
of line widths comparable to the range of black hole masses. 

Alternatively, and more controversially, it may be that 
the BLR is not virialised at all and the virial relations only 
prove accurate due to their luminosity dependence. The 
SMBH mass of a QSO and its luminosity must be strongly 
covariant and the ability of virial estimators to determine 
Mbh rnay be due to this simple relation. The line width term 
may be redundant. Indeed we (and other authors e.g. Cor- 
bett et al. 2003; Shen et al 2007) find the average line width 
changes very little over the luminosity range of our sample. 
To be consistent with the virial mass estimators this requires 
that both SMBH mass and Eddington ratio for QSOs vary 
as luminosity to the power of ~ 1/2. A conspiracy which 
must be viewed with some caution. 

The reason for the lack of observed scatter in broad line 
widths is uncertain, but it is apparent that this observation 
has consequences for black hole mass estimation in AGN. 
While this may imply that virial masses are not necessarily 
unbiased estimators of Mbh there is evidence from other 
studies that they can be used as an indicator of the mass of 
a black hole. 
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Figure 13. A sketch showing the assumed geometry of an AGN. 
The central SMBH and BLR are surrounded by an obscuring 
torus which constrains the allowed observing angle (9) to be less 
than some opening angle {dm) 



6.4 Constraining BLR geometry 

The low dispersion in broad-line width found for luminous 
QSOs is also a strong constraint on the velocity field of the 
BLR. 

In all but the case of a spherically symmetric BLR the 
velocity dispersion we measure from a spectrum will depend 
on the viewing angle to the QSO. Given a model for the 
BLR we can calculate the expected dispersion in measured 
line widths due to variations in observing angle. If the ex- 
pected dispersion is greater than we find in our data we can 
rule out that model for the BLR. Since we take no account 
of the myriad of other sources of scatter in the line width 
distribution, this equates to a very strong constraint on a 
given BLR model. 



6.4-1 Pure planar/polar BLRs 

The first models we consider are planar velocity fields in 
which all the velocities are confined to the z = Q plane 
(cylindrical coordinates), and polar fields in which all the 
velocities are in the z direction. In these cases we expect 
our measured IPV widths to be modified by 
Planar: log(/Py) oc log(sin(6l)) 
Polar; log(/PF) oc log(cos(6l)) 
where 9 is our observing angle to the QSO (see Fig. 13). 
Assuming an opening angle (Sm), which we take to be con- 
strained by an obscuring torus (Fig. 13), we can calculate 
what the expected dispersion in measured line widths would 
be due solely to orientation effects for planar and polar ve- 
locity fields. 

Fig. 14 shows how much scatter in line widths we expect 
in each case. On both plots we indicate with a dashed line 
0.07 dex dispersion which is the lowest we measure in our 
sample (Fig. 7). 

In the planar case the expected dispersion is ~ 0.2 dex 
for all opening angles. We find somewhat less dispersion in 
IPV widths than this in all but two of our L — z bins. Hence 
the velocity field the BLR cannot solely be planar but must 
include some other component. 

For the polar velocity field the expected dispersion is a 



strong increasing function of 6m- For modest opening angles 
the resulting dispersion is low and it is not until 6m ~ 55 ° 
that this becomes comparable to the dispersions we measure 
in our data. Hence if the BLR is characterised by polar flows, 
the opening angle to luminous quasars must be less than this 
value. 



6-4-2 Planar/Polar BLRs with a random/ spherically 
symmetric component 

It is not likely that in any model for the BLR all the 
velocities are confined to a single plane/direction. Instead 
there will always be some random component to the veloc- 
ity field and we include this in our models with the following 
parametrisation: 

Planar: log(/PV) oc log(a sin(6l) + (1 - a)) 
Polar: log(JPV) oc log(a cos(6i) + (1 - a)) 

In these models the trigonometric term represents a ge- 
ometrically constrained component to the BLR velocity, and 
the (1 — a) term represents a random/spherically symmetric 
component. Hence a model with a = 1 represents the pure 
planar/polar cases outlined above, and if a = the BLR 
velocity field is spherically symmetric. In these cases, given 
values for a and 6m we can again calculate the expected 
dispersion in log(/PV) to compare with our results. These 
calculations are illustrated in Fig. 15. 

In this figure the grey scale indicates the expected scat- 
ter at a point in the a — 6m plane. The contours are at 
increments of 0.035 dex, and so the second contour (heavy 
line) represents the 0.07 dex scatter we observe for the most 
luminous QSOs in our sample. Hence, for these luminous 
objects the region of parameter space above the solid line is 
ruled out in our analysis. 

In the planar case, for which the constraints are 
stronger, taking a believable value for the opening angle 
(~ 45°) our data show that roughly half of the contribution 
to the velocity field of the BLR must come from a symmetric 
component. 



6-4-3 A hybrid BLR 

As a final model we consider a BLR with both planar and 
polar components to the velocity field as well as a spherically 
symmetric component. We model this with: 

log(JPV) oc log(asin(e) + hcos{6) + {l-a-h)) 
And so for a given (a, 6, dm) we can calculate the expected 
dispersion in measured line width. Fig. 16 shows the expected 
dispersion in the a — b plane for differing opening angles. In 
each plot the origin (a = & = 0) represents a spherically 
symmetric BLR and the line a + b = 1 represents a BLR 
with no symmetric component. As one moves from left to 
right in each plot the BLR becomes more planar, and from 
bottom to top the velocity field is more polar. 

As expected, for small opening angles our data are not 
strong constraints on the velocity field of the BLR, and only 
the most disk-like models are rejected. For larger opening 
angles there is a larger region of the parameter space ruled 
out by our data, predominantly the very planar or very polar 
models. 
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Figure 14. These plots show how the dispersion in measured hne width is dependent on the opening angle of the source in the special cases 
where either all of the BLR velocities are in the planar or polar directions. On each plot the dashed line shows the dispersion = 0.07 dex 
level, the lowest dispersion in line width we measure in our data. 
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Figure 15. Here we demonstrate how adding a spherically symmetric component to the BLR velocity field affects the expected dispersion 
in broad line widths. Again we assume that all of the other velocities are either in the planar or polar directions. In our parameterisations 
a ~ indicates a BLR dominated by a spherically symmetric velocity field, while a ~ 1 indicates no symmetric component equivalent to 
the models shown in Fig. 14. The grey scale and contours indicate the expected dispersion. Contour levels are in increments of 0.035 dex 
and the second contour (heavy) shows the 0.07 dex scatter we measure for the most luminous objects in our sample. The parameter space 
above this is inconsistent with our analysis. 



7 CONCLUSIONS 

We have measured the 50 % IPV width of the Mg ll line 
in QSO spectra from the SDSS, 2QZ and 2SLAQ surveys 
and find a strong correlation between the dispersion in IPV 
widths and the optical luminosity of QSOs. If we assume 
there exists a virial relation of the form of equation 2 this 
implies that there is an equivalent reduction in the disper- 
sion in Mbh- On face value this is in excellent agreement 
with models for the QSO population proposed by Hopkins 
et al. (2005b). 

However, the remarkably low scatter we find in our IPV 
width measurements, in particular for the more luminous 
objects, has implications as to the validity of virial mass 
estimators. We find less scatter in IPV widths of luminous 
QSOs than is intrinsic to the r — L relation. While the r — 



L relation is not defined for these very luminous objects 
these results are at odds with the practise of virial black 
hole estimation. Possible explanations for this observation 
include very low scatter in the r — L relation for bright QSOs 
or a BLR which is not fully virialised. In either case it is clear 
we are yet to gain a full understanding of the process of virial 
black hole mass estimation and one must be cautious when 
performing or interpreting these estimates. 

Finally we show that the observed scatter in IPV widths 
can be used to constrain models for the velocity field of the 
BLR. We show how variations in observing angle to sources 
affects the scatter in measured IPV widths for a series of 
simplified geometries and derive constraints on these models 
from our data. Perhaps of most interest is the rejection of 
a pure planar BLR regardless of the assumed opening angle 
to the source. 
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Figure 16. These plots show the expected dispersion in measured line widths for our composite BLR model which include planar, polar 
and symmetric components. Each plot shows the expected dispersion (grey scale + contours as in Fig. 15) for a given AGN opening angle 
{9m)- In our parametrisation as a increases the BLR is more planar, and as b increases more polar. (0,0) indicates a totally symmetric 
velocity field while all models on the line a + b = 1 have no symmetric component. 



A table containing the Mg ll IPV widths for the quasars 
studied in this work is available through the 2SLAQ website 
(WWW.2SLAQ.INF0). 
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